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ABSTRACT | 


This program analyzes the errors of radiolocation predictions and evaluates 
the average range errors, the root mean square range error, and ten, twenty, and 
thirty minute cumulative averages of the range errors for both the following: 

a) Range errors vf the predicted range with respect to the true range 
of the transmitter where range is measured from the receiver. 
b) Range errors of the predicted transmitter position where range is 
measured from the known transmitter location. 
In addition to this, the computer punches up paper tapes for operating an X-Y 
plotter. Twotypes of graphs are obtained in this manner. One is a plot of 
predicted range error vs. time where range is measured fromthe receiver. The 
other is a plot of the range and bearing of the predicted transmitter position 
with respect to the known transmitter location. The latter graph is referred to 
as a target plot. This program is coded inSPAR. The computer involved in 


this case is a Control Data Corporation G-20A using the SPACE programming 


system. 


1. PROBLEM STATEMENT 


Measurements of the ver‘ical angle of incidence and azimuth of an arriving 
HF radio signal are made with an interferometer. The vertical angle of 
incidence and vertical incidence 1onogram data are used to predict the range 
Or range error in the case where the location of the transmitter 1s actually 
known. The azimuth, time, and interferometer predictions from the ionogram 
studies and other information such as the propagation mode, provide the input 
to this program. The program analyzes the errors of radiolocation predictions 
and evaluates the range errors, the root mean square range error, and ten, 
twenty, and thirty minute cummulative averages of the range errors for both 
the following: 


a) Range errors of the predicted range with respect to the true range 
of the transmitter where range is measured from the receiver. 


b) Range errors of the predicted transmitter position where range is 
measured from the known transmitter location. 


In addition to this, the computer punches up paper tapes for operating an X-Y 
plotter. Twotypes of graphs are possible. The first is a plot of the range 
errors mentioned in (a) above as a function of time (see Hgure 3). The second 
is a plot of the range and bearing of the predicted transmitter position from the 
true transmitter location (see Figuse 6). This is referred to as a target plot. 
In addition calculations are made to obtain the range errors described in (a) and 
(b) above which are then averaged according to the following scheme: 


a) Ten minute cumulative averages of the range error starting on every 
ten minutes beginning with the period during which data was first 
taken. 


b) Same as a) except for twenty minutes. 
c) Same as a) except for thirty minutes. 
d) Average range error. 


e), Root-mean-square range error. 


‘ 


As far as the 10-, 20-, and 30-minute cumulative averages are concerned, an 


hour is divided into fixed time intervals: i.e., 6-10 minute intervals, 3-20 


minute intervals, and 2-30 minute intervals. These intervals are built into the 
program and are not available to the operator. For example, consider that the 
first data was recorded at hour X+ 45 minutes. The printout starts with hour 
X. The 10 minute cumulative average starts with X + 40 as does the 20 minute 
cumulative average. The 30 minute cumulative average starts with X + 30. 
As was previously implied, the cumulative averages and final averages are printed 
by a line printer. Since the graphs are plotted on blank paper, the proper 
coordinates must be plotted on the paper also. This 1s done by separate paper 
tapes whose use is described later under Operating Instructions. Because the 
coordinates are plotted according to a certain scale, it is necessary that the 
data also conform to this scale which is as follows: 


a) Range Error vs. Time Plots: 50 kilometers and one hour are 
equivalent to 6666 and 8000 plotter units respectively. 


b) Target Plots: 100 kilometers are equivalent to 8888 plotter input 
units. 


The target plot has it origin at the center of the graph paper, therefore, the 
plotter is zeroed at this point (see Figures 4 and 5). It is not necessary to 
change the calibration or zero ppint of the X-Y plotter for the Range Error vs. 
time plot as the horizontal centerline of the graph represents 0 kilometer range 
error and the vertical centerline represents the beginning of the second hour as 
the time axis of the Range Error plot covers two hours. The calibration for 
the X-Y plotter is discussed in the operating instructions for this program. A 
total day's data may easily run more than two hours in which case the data is 
plotted on succeeding sheets of graph paper. There is, however, no overlapping 
of the graphs. Note that the odd hours have negative X components but this is 
resolved by the program. Also inthis general problem statement, it is worth 
while to point out some of the contingencies of the data which have been designed 
into the program. These are as follows: 


a. Target Plot: Sometimes the point to be plotted exceeds some 
maximum range (greater than + 75 kilometers north and south and 
greater than + 100 kilometers east and west) with respect to the 
true transmitter position. Inthis case, this particular point's 
range and bearing plus some related information is printed out 
and not plotted. This point, however, is included in the cumulative 


averages. 


b. Range Error vs. Time Plot: Sometimes the range error exceeds 
+ 50 kilometers in which case the error aid time are printed along 
with some related information and is not piotted. However, it is 
included in the cumulative averages. 


c. Raw Data: Sometimes there are no range predictions available when 
this data is required. In this event, the data card is read but ignored 
on input. The reason that data cards of this nature are in the data 
decks is that part of the i::put data will be punched up and used by 
programs, other than this one, which are still in the developmental 
stages. 


d. An identification, which labels the data following it, is punched in 
the paper tapes prepared by the computer. The tapes are punched 
with each class following sequentially and in the following format. 
For each class: eighteen inches of leader, the data class identifi- 
cation, six inches of blank tape, and then the data block itself 
followed by a stop code for the plotter. 


e. If no data exists for a particular class or graph (i.e., Range error 
data available for lst and 2nd hours, none for third and fourth, but 
finally some for the fifth hour) appropriate messages are punched 
On paper tape. 


f. Printout, Cumulative Averages: The value of the average, the number 
of datum used in computing the average, and the time of the last datum 
used in computing that average are printed for the ten, twenty, and 
thirty minute averages in a suitable format. 


2. COMPUTATIONAL METHODS 


With the exception of the cumulative averages no equations are discussed in 
this section. The bulk of the arithmetic is performed by the RANGE and COORD 
subroutines described by RRL Publication No. 231. The remaining mathematics 
consist of simple trigonometric functions which are readily apparent from the 
flowchart and listing of this program. Therefore, most of the program involves 
manipulation of the data according to certain techniques which are outlined and 
described in this section. Basically the raw data 1s read in and stored on magnetic 
tape which is then scanned seven times for each type of plot. The data is scanned 
once for computing the cumulative average and then once for each of the six data 
classes which are determined by the antenna pair and propagation mode associated 
with each sub-group of raw data. 

The program is divided up into five sections whose functions are described 


as follows: 


A. Section I 


I. 


2. 


Reads two cards of initial data whose contents are described later 
under Input and Output. 


Converts the longitude and latitude of the receiver and transmitter 
from degrees, minutes, and seconds into degrees and fractional 
degrees which is then converted into radians. 


Computes the true range between the receiver and the transmitter 
which is stored for future reference using the RANGE subroutine. 


Determines the parity (even or odd hours) of the reference time of 
the initial data. This reference time becomes the basis of a relative 
time reference hour. The time in hours and fractional hours which 
appears in the various program printouts in time relative to the 
initial time reference, .i.e., time reference of 2100 hours anda 
relative time of 1.75 hours imply that the true time is 21:45:00; a 
time reference of 2200 hours and a relative time of 3.75 hours imply 
that the true time is 23:45:00 


Reads in the raw data, packs this in a 256 word temporary holding 
zone, which when full is dumped onto magnetic tape. When the last 
of the raw data is read, an "end'' code is placed behind the final 
data all of which is stored on magnetic tape. As the data is read in, 
the following tests and conversions are made: 


a. If the predicted range is to be used and is zero, this data 
card is thrown out. 


b. Time is converted to fractional hours and added to the either 
odd or even time reference and stored thus. At each hour 
change, the time reference is incremented accordingly. 


c. If time skips an hour, this fact is processed and the time 
reference is incremented accordingly. 


B. Section Il 


1. 


Prints the program heading and the proper heading for the cumulative 
averages currently being computed. 


Checks to see which plot type is currently being analyzed and on the 
basis of this proceeds to either (a) or (b) as follows: 


a. Target Plot: computes the distance between the predicted 
transmitter position which is then fed, along with the current 
time, to the cumulative average machinery. The bearing of 
this range is computed here also but is not used. To obtain 
these results the program proceeds as follows: 


1. If range error is given, it is added to the true range 
between the receiver and the transmitter to obtain the 
predicted range. 
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2. If the predicted range between the receiver and trans- 
muiiter 1s given, it is used directly. 

3. The COORD subroutine 1s now used to obtain the 
coordinates of the predicted transmitter pointer location. 

4. These coordinates alung with the coordinates of the true 
transmitter position are fed into the RANGE subroutine, 
the output of whichis the desired distance between the 
predicied and true transmitter locations 

5. The coordinates of the true transmitter location are then 
subtracted from the predicted transmitter position and 
the relative longitude and latitude obtained are used to 
obtain the bearing of the range Computed in step No 4 
(See program flowchart and Figure | tur details) 


b. Range Error Plot determines the range error ahich. along 
with the current time, is fed to the cumulative average machinery. 
If the range error is given nothing needs to be done, the range 
error is used directly If the predicted range 15 given, the 
predicted range minus the true range gives the desired range 
error. (See program flowchart and Figure 1 for details) 


Scans the raw data on magnetic tape. computes the Cumulative averages 
and prints the results ina suitable format as they are computed The 
cumulative averages are determined by 


aa , - 
ang Xx where X is the cumulative averape (') 
ee ae 1 : 
i-l 
A ten-minute average involves all the points within a ten minute intertad 


Tren it 1s reset and started over The same applies to the twenty and 


thirty-minutes averages. 


Computes a total average and RMS average which are printed out along 
with the total number of points used in computing these averaypes. The 
total average 1s computed overthe whole data run trom (1) ands not 
reset at the end of every specified time interval. The tinal RMS averape 
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This final equation is the one solved by the program. 


Section III 


1. Punches out the required amount of leader and the proper identification 
for each class of data. 

2. Contains the machinery for scanning the raw data picking up only the 
data required by the current class. If none is found, the message 
"NO DATA OF THIS TYPE" is punched on the paper tape along with 
a stop code for the plotter. 


3. Transfers to the next section, if good usable data 1s picked up in the 
raw data scan, where it is processed. 

4. Punches remainder of data contained in a punch zone (described later) 
plus a stop code for the X-Y plotter when the end of the raw data is 
detected. 


5. Sets up the proper paper tape identification for the next data class to 
which it then proceeds to scan. 


6. Transfers to an exit zone (See Section V), after scanning the last data 
class, which determines whether another plot is desired. If so, the 
program is restarted through channel two (Range error plot) and if 
not, terminates the program. 


Section IV - checks to see which plot type is currently being processed 
and on the basis of this proceeds to either (1) or (2) as follows: 


l. Target Plot 


a. Computes the values of range and bearing as des« ribed in Section 


II paragraph 2a. 


b. Checks the values of the X and Y components of the range and bearing 
of the predicted transmitter location to determine whether or not 
the X and Y components are within certain limits. If these limits 
are exceeded, the data is printed out and not processed beyond this 
point. 


Scales the proper values of the X and Y components, assuming that 
they fulfilled the conditions of the previous step, forthe X- Y 
plotter and restores these numbers into locations common with the 
Range error channel of this section and exits to the next section. 


Cc. 


2. Range Error Plot 


a. Computes the value of the range error as described in Section ITI, 
paragraph 2b. 


b. Checks the value of the range error to determine whether or not 
the range error exceeds a certain limit. If so. the data is re- 
jected and if not. the program proceeds toc. 


c. Determines upon which two-hour- graph the first time for a particu- 
lar graph falls. For each two-hour-graph that there is no data to 
plot, the message ''NO DATA THIS GRAPH!" is punched into the 
paper tape followed by a stop code for the X - Y plotter. 


d. Determines the parity of the relative time in order to establish 
an odd hour time reference and an even hour time reference. 
This step and the preceding one (c & d) are executed only for the 
first time through on each particular data type. 


e. Compares the current relative time with the odd and even time 
references. If the current time is odd, it is converted intoa 
suitable negative ''time" which fulfills previously mentioned 
criteria. If the current time is even. no conversion is made. 
This section also contains the machinery for incrementing the 
odd and even reference time as necessary. Machinery is also 
built in to determine if the current time skips a two-hour~-graph. 
If this is true, a suitable message (‘'NO DATA THIS GRAPH") is 
punched into the paper tape followed by a stop code. (Details are 
evident from the flowchart and listing. ) 


f. Scales the proper values of range-error and time for the X - Y 
plotter and restores these numbers into locations common with 
the Target Plot channel of this section andexits to the next 
section. 


Section V 


1. Receives the number pairs from Section IV and breaks these down intoa 
a sign character and single integers stored in sequential memory locations. 
(As a result of the scaling done in section IV, the numbers coming into 
Section V are in the form tdddd, a sign and four digit integer. ) 


2. Converts the integers into the corresponding characters required by the 
X - Y plotter. 


3. Formats these characters for the X - Y plotter. 


4. Contains the machinery for accumulating the character strings in a punch 
zone which when full is dumped onto paper tape. 


5. Determines finally whether or not a second plot is desired and then either 
exits to SPACE or restarts the program with the running averages 
(Section II). 
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In order to operate the X-Y plotter to piot one point, the following information 


is needed and is entered in the order snown. 


1. X follows 
2. Sign 
Th igi 
ousands digit X coordinate 
Hundreds digit 
Tens digit 


Units digit 


Y follows 
sign | 


Thousands digit Y coordinate 


o> EN °° © Sn © DE ©. - 


Hundreds digit 

ll. Tens digit 

12. Units digit i 

13. Plot instruction instructs plotter to plot this point. 


Essentially, then, Section V takes two four-digit 1utegers plus signs and translates 
these into the above format. In the event that one of the digits «an not be translated 
into the corresponding plotter character, the message ' TRANSLATOR ERROR" 


will appear on the line printer and that data 1s rejected. 


3. LIMITATIONS AND ACCURACY 


Although the arithmetic of the two subroutines. RANGE and COORD, 1s double 
precision, the arithmetic of this control program 1s all single precision which 
limits the accuracy to six decimal digits plus a two digital exponent. This 
accuracy is more than sufficient since the resolution of the plotter 1s limited 
to four significant digits and most of the input data has the same accuracy. Other 
limitations are discussed inthe report on the subroutines, RRL Publication No. 
231, but normally these will not be of interest. This program assumes the tol- 


lowing sign conventions on longitude and latitude: 


a. longitude east of Greenwich is positive 

b. longitude west of Greenwich is negative 

c. latitude north of the Equat>::s positive 

d. latitude south of the Equatoris negative 
In the operation of the X-Y plotter using the proper tapes prepared by this program, 
there is no stop code after the second hour of each two hour plot of the Range 
Error vs. Time plot. Therefore, when this hour change occurs the operator must 
stop the plotter's paper tape reader and change the sheet of graph paper. other- 
wise the third and fourth hour, etc. will be plotted over the first and second 
hour's data, etc. This program will process a data run of up to eight hours 


duration; this is more than anticipated at the present. 


4. COMMON STORAGE AND INDEX REGISTERS 
This program requires fourteen memory locations of common storage which 
are used for communicating with the RANGE and COORD subroutines. The program 


uses five index registers but no common index registers. 


5. OPERATING INSTRUCTIONS - 1/O DESCRIPTIONS 
The program reads in two cards of initial data which contain the following 


information. 
A. First card 


1. Degrees 


2. Minutes Receiver longitude 

3. Seconds 

4. Degrees 

5. Minutes Receiver latitude 

6. Seconds 

7. Time Reference - Hour in which data begins 

8. Time Reference - Some time in minutes equivalent to one minute 


less the time on the first data card (i.e. if time on first card is 
5.06 minutes, this time reference should be 4.06, etc.) 


9. Plot - type of plots desired from the program: 
This will be an integer 1, 2, or 3 signifying: 


Target plot only 
Range Error plot only 
= both Target plot and Range Error plot in that order. 


1 
2 
3 


B. Second card ‘ 


lL. 
é. 
3. 
4. 
5 

6. 


Degrees 

Minutes Transmitter longitude 
Seconds 

Degrees 

Minutes Transmitter latitude 


Seconds 


The remaining cards, which vary in number, contain the raw data which 1s to 


be processed. 


Azimuth a 

Time (in minutes) Punched up by computer interferometer program. 
Control Number A 

Predicted Range 

Range Error 

Control Number B 

Control Number C punshal br huuan operate 
Letters SKIP if time s'ps an hour 


Letters END to denote end of raw data 


The program uses either Prec:cted Range or Range Error. if Predicted Range is 


used, enter 0.0 in range error field; if Range Error is used, enter 0.0 in 


predicted range field. 


The Control Numbers have the following significance: 


A 


oO www PP 


1: 0 - 120 Interferometer Antenna pair type data. 
2: 0 - 240 Interferometer Antenna pair type data. 
1: Propagation mode one hop E type data. 

2: Propagation mode two hop E type data. 

3: Propagation mode one hop F type data. 

1: Predicted Range to be used for this data. 


2: Range Error to be used for this data. 


1) 


Normally SKIP and END are absent from these cards. SKIP 1s punched in 
columns 73-76 of the data card only if time skips one hour. END punched in 
columns 73-80 denote the end of the raw data. Upon detecting END on a data card. 
the computer will begin to process the data read in prior to this card. See Figures 
2a and 2b for data card formats, 

The calibration procedures for the X Y plotter are adequately explained in 
two instruction manuals supplied with this instrument. For the purposes of this 
program, the plotter is zeroed in the center of a sheet of graph paper and the 
maximum length available for the positive ordinate axis 1s made equivalent to 
6666 plotter units which implhes that the maximum length of the negative ordinate 
axis is -6666 plotter units. The length of the ordinate axis 1s then made equivalent 
to 6666 plotter units on the abscissa (See diagram below). This completes the 
calibration of the X-Y plotter for both the Target and Range Error vs Time plots. 
It may be necessary to check this calibration from time tu time during the course 


of plotting a large number of graphs 


CALIBRATION DIAGRAM 


11 x 1? praph 
paper 


bobb 


The coordinates are then plotted using an appropriate paper tape followed by 


the actual data. 


Where: K 


Fipure l, 
T 
dl 


PT 


'Z 
RE 


The predicted 


ls the receive. position, 

is the transmitter position, 

rs the azimuth of the incoming radiowave, 

is the predicted transmitter position, 

is the distance between the predicted and the 


true transmitter position, 


is the bearing of E, 
is the range error between the true and the 
predicted transmitter position, 

range and true range are as noted above, 
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Figure 2a, Initial Data Cards (only two required for one 
data run,) 
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